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摘 要 : 土地 沙漠 化 是 当今 全 球 最 严峻 的 环境 问题 之 一 ,也 是 我 国 西北 干旱 区 
境 问题 。 开 展 土地 沙漠 化 动态 研究 有 助 于 沙漠 化 灾害 防治 与 国土 资源 空间 
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临 的 主要 生态 环 
发 。 以 2000 一 2018 


年 MODIS-NDVI1 km 月 数据 为 基础 ,运用 栅 格 累加 法 与 转移 矩阵 法 对 塔里木 盆地 土地 沙漠 化 的 变 
化 速率 、 空 间 特 征 及 其 成 因 进 行 探索 分 析 , 结 果 表 明 :(1) 土地 沙漠 化 是 一 个 具有 年 际 波动 和 累积 
效应 的 动态 变化 过 程 ;(2) 栅 格 累加 法 具有 较 强 的 连续 动态 分 析 能 力 ,能 准确 分 析 土 地 沙漠 化 变化 
速率 演变 趋势 .变化 的 空间 特征 ,可 以 定量 分 析 土 地 沙漠 化 原因 ;(3) 塔里木 盆地 沙漠 化 退化 主要 
是 草地 的 退化 ,沙漠 化 改善 则 主要 是 林地 和 耕地 的 改善 ;(4) 2000 一 2018 年 ,塔里木 盆地 土地 沙漠 
化 整体 呈 极 重度 沙漠 化 逐 级 向 非 沙漠 化 变化 的 改善 趋势 ,转移 路 径 为 : 极 重 度 沙 漠 化 一 重度 沙漠 化 
一 中 度 沙 漠 化 一 轻 度 沙漠 化 一 非 沙 汉化 ,改善 区 主要 分 布 在 绿洲 区 ,恶化 区 主要 在 塔里木 河中 游 
与 车 尔 臣 河 下 游 ,土地 沙漠 化 成 因 主要 为 人 类 活动 。 此 文 研究 方法 为 土地 沙漠 化 动态 监测 提供 了 


一 种 新 思路 ,其 研究 结果 可 为 塔里木 盆地 国土 资源 


开发 和 沙漠 化 治理 提供 理论 依据 。 
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文章 编号 : 


土地 沙漠 化 是 当今 世界 最 严重 的 环境 问题 之 
一 ,影响 全 球 2/3 的 国家 和 10 余 亿 人 口 "。 我 国 受 
土地 沙漠 化 影响 严重 ,沙漠 化 面积 占 国 土 面积 达 
27% ”。 新 疆 塔 里 木 盆地 尤为 严重 ,每 年 沙漠 向 绿 
洲 推移 5 m 以 上 ,1296 km 土地 遭受 沙化 威胁 ” , 农 
田 、 公 路 .铁路 等 遭受 沙 埋 ,风沙 灾害 频 发 ,严重 影 
响 塔 里 木 盆 地 甚至 我 国 西北 地 区 的 人 民生 产 和 生 
活 ,路 需 对 塔里木 盆地 土地 沙漠 化 动态 变化 开展 
研究 。 

植被 是 联系 土壤 、 水 体 和 大 气 之 间 的 关键 纽 
带 ,是 土地 沙漠 化 动态 变化 监测 的 重要 指标 “ ,是 遥 
感 监测 土地 沙漠 化 的 “指示 融 ””。 归 一 化 植被 指数 
(Normalized difference vegetation index, NDVI) 能 较 
HE J We ie Te He oe BE , LE Fie hs DN Ps Az 
内 外 学 者 运用 Mann-Kendall 检 验 法 ,趋势 分 析 法 、 
主 成 分 分 析 法 、 差 值 法 等 方法 对 NDVI 动 态 变 化 进 
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行 了 大 量 监测 研究 ™” ,其 中 差 值 法 计算 简单 结 
直观 被 广泛 采用 一 六 。 然 而 ,大 部 分 研究 仅 用 几 期 
遥感 影像 进行 差 值 分 析 ,其 监测 结果 可 能 与 实际 情 
况 存在 一 定 差 异 。 植 被 盖 度 受到 降水 量 ` 气温、 人 
类 活动 等 因素 影响 , 仅 用 几 期 影像 进行 监测 存在 偶 
然 性 偏差 ,更 换 不 同年 份 的 影像 所 得 结果 可 能 不 相 
同 。 栅 格 累加 法 以 长 时 间 序 列 数据 为 基础 进行 差 
值 计算 可 以 避免 此 类 偶然 性 缺陷 ,能 真实 地 反映 
植被 变化 。 

本 文选 择 2000 一 2018 年 月 尺度 的 MODIS ND- 
VI 数据 (MOD13A3), 采 用 栅 格 累 加 法 和 转移 矩阵 
法 对 塔里木 盆地 土地 沙漠 化 变化 进行 研究 ,以 期 得 
到 塔里木 盆地 准确 的 土地 沙漠 化 变化 过 程 与 特 
征 。 其 研究 结果 为 塔里木 盆地 国土 资源 的 开发 利 
用 做 数据 支撑 ,为 “丝绸 之 路 经 济 带 ” 相 似 区 域 提 供 
参考 。 
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1 研究 区 概况 


塔里木 盆地 是 典型 的 内 陆 盆地 ,盆地 降水 稀 
> 洪 在 蒸发 大 ,植被 稀 玻 ,属于 极端 干旱 区 。 盆 地 
年 均 降 水 量 <50 mm ,潜在 蒸发 量 达 3400 mm ,降水 
难以 形成 有 效 地 表 径 流 ,难以 满足 植被 正常 生长 需 
求 , 植 被 主要 依靠 河流 水 和 地 下 水 生长 。 盆 地 主要 
河流 有 塔里木 河 .阿克苏 河和 了 田 河 、 开 都 河 、 喀 什 
噶 尔 河 .车 尔 臣 河 (图 1)。 自 然 植 被 大 多 沿 河 分 布 ， 
ERA HH (Populus euphratica ) EH (Tamarix chi- 
nensis ) IEW (Alhagi sparsifolia ) EE% ( Karelin- 
ia caspia ) 芦苇 (Phragmites communis)", Ab 
ERKE HKW E , BAA VURE Z ,质地 较 粗 ， 
易 被 风蚀 。 


0 150 km 


图 1 塔里木 盆地 示意 图 
Fig. 1 The sketch map of Tarim Basin 


盆地 人 口 1095.96x10', 人 类 活动 较为 活跃 。 随 
着 人 口 的 增多 ,耕地 面积 逐步 扩大 ,河水 与 地 下 水 
过 度 引 用 ,河道 断 流 和 地 下 水 位 降低 ,下 游 河道 廊 
道 植 被 退化 ,导致 沙漠 化 发 展 迅速 。 不 科学 灌溉 耕 
地 形成 士 壤 盐 渍 化 ,大 面积 耕地 被 迫 扬 荒 ,进一步 
加 剧 土地 沙漠 化 。 虽 然 塔 里 木 河 输 水 工程 取得 了 
巨大 的 成 效 , 但 塔里木 倪 地 的 沙漠 化 依然 严重 "1。 


2 数据 来 源 


MODIS-NDVI 数 据 质量 稳定 .周期 性 强 、 容 易 获 
Be) Be YZ AA), AS SC HH ae H Se EE] USGS 
Ke AG HS MODIS 植被 产品 (MOD13A3 ) ,空间 分 辩 率 1 
km。 本 文 根 据 研究 区 需要 ,选择 轨道 号 为 h23v04、 


h23v05, h24v04, h24v05. h25v04, h25v05, 2000— 
2018 年 月 合成 数据 ,共计 1362 景 。 为 进一步 分 析 土 
地 沙漠 化 成 因 , 从 http:/www.dsac.cn/ 下 载 新 疆 2010 
年 土地 利用 图 ,空间 分 辨 率 为 1 km。 


3 研究 方法 


3.1 植被 盖 度 计算 

运用 MRT(MODIS reprojection tool ) 对 影像 进行 
拼接 .裁剪 .投影 转换 ,得 到 塔里木 盆地 2000 一 2018 
年 NDVI 月 数据 集 , 值 域 为 [-3000, 10000]。 为 消除 
KAS .太阳 高 度 角 、 季 节 等 因素 干扰 ,对 NDVI 
月 数据 进行 年 最 大 值 合成 ,生成 年 最 大 值 NDVI 

NDVI, = max(NDVI,), j=1,2,3,.…,12 (1) 
式 中 :NDVL 代 表 第 i 年 的 NDVI 最 大 化 合成 数据 ; 
NDVJ; 表 示 第 i 年 ) 月 的 NDVI 数 据 。 

影像 中 NDVI 值 小 于 -1000 的 区 域 为 水 域 ,为 消 
除 水 域 变化 对 沙漠 化 变化 的 影响 ,将 其 小 于 -1000 
的 NDVI 栅 格 值 赋 为 10000。 影 像 中 0 到 -1000 的 区 
域 为 无 植被 区 ,NDVI 赋 值 为 0。 

NDVI 是 利用 植被 光谱 特征 计算 出 来 的 一 种 
植被 指数 ,缺乏 明确 的 物理 意义 ,无 法 进行 地 面 验 
证 "。 本 文采 用 线性 像 元 二 分 模型 将 NDVI 转换 成 
植被 覆盖 度 ,计算 公式 为 : 

_ NDVI-NDVI,,, 

“ NDVI, -NDVI a 
IP FW HORE E BE; NDVT。 为 植被 盖 度 接近 
100% 的 NDVI 值 ;NDVLa 是 无 植被 地 的 NDVI 值 。 
NDVL. 选 择 每 年 NDVI 值 都 接近 接近 10000 的 区 域 ; 
NDVILs 选 择 无 植被 的 塔克拉玛干 沙漠 流动 沙丘 。 
每 年 选择 相同 NDVI 与 NDVI 区 域 减少 土壤 背景 
值 对 计算 精度 的 影响 。 为 方便 计算 ,将 像 元 值 域 线 
性 调整 到 [0, 100] ,得 到 2000 一 2018 年 的 植被 盖 度 
年 度数 据 集 。 
3.2 土地 沙漠 化 分 级 

本 文 参考 王涛 提出 的 沙漠 化 分 级 体系 局 ,结合 
塔里木 盆地 实际 情况 进行 分 级 ,将 研究 区 土地 沙漠 
化 分 为 非 沙漠 化 、 轻 度 沙 漠 化 .中 度 沙 漠 化 .重度 沙 
漠 化 、 极 重度 沙漠 化 5 级 ( 表 1)。 

3.3 栅 格 累加 计算 
土地 沙漠 化 受 自然 、 人 为 等 因素 影响 而 呈 波 动 


(2) 
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表 1 土地 沙漠 化 分 级 


Tab. 1 Desertification classification in Tarim Basin 


沙漠 化 分 级 非 沙 漠 化 轻 度 沙漠 化 中 度 沙漠 化 重度 沙漠 化 极 重 度 沙漠 化 
植被 盖 度 /% [100.0,60.0) [60.0 ,40.0) [40.0,26.7) [26.7,13.3) [13.3, 0) 


变化 ,并 具有 一 定 累 加 效应 ,本 文 运用 栅 格 累加 法 
对 多 年 土地 沙漠 化 进行 研究 ,分 析 多 年 土地 沙漠 化 
变化 状态 以 及 空间 格局 ,计算 公式 如 下 : 

E,= YW, -U,)/n (3) 


式 中 :为 第 n 年 土地 沙漠 化 累加 效应 ;0 为 第 n 年 
土地 沙漠 化 现状 ;00 为 基准 年 土地 沙漠 化 现状 。 

栅 格 累加 的 实质 是 计算 多 年 各 个 栅 格 像 元 植 
被 盖 度 的 年 平均 变化 速率 , 正 值 表 示 栅 格 像 元 植被 
盖 度 增加 , 负 值 表示 减少 ,数值 的 大 小 表示 平均 变 
化 的 快慢 。 依 据 土地 沙漠 化 的 分 级 标准 实地 调查 
数据 分 析 和 反复 尝试 ,最 终 栅 格 累加 分 级 标准 如 表 
2 所 示 。 
3.4 土地 沙漠 化 成 因 分 析 

计算 相 邻 2 期 土地 沙漠 化 等 级 图 进行 转移 算 
阵 , 分 析 不 同 沙 漠 化 等 级 之 间 的 转移 量 与 空间 分 
布 ,结合 土地 利用 类 型 分 析 土 地 沙漠 化 成 因 。 前 人 
研究 显示 ,塔里木 盆地 土地 沙漠 化 主要 由 水 资源 不 
Fr BA FA a E eae CO ,过 度 放 牧 、 樵 柴 毁 林 ALAS 
中 草药 相对 较 少 。 因 此 在 分 析 土 地 沙漠 化 变化 成 
时 ,将 耕地 建设 用 地 上 的 土地 沙漠 化 变化 划 为 
人 类 活动 直接 因素 ,草地 林地 水域. 未 利用 地 上 
的 土地 沙漠 化 变化 划 为 人 类 活动 间接 因素 。 根 据 
转移 矩阵 计算 结果 ,统计 分 析 土 地 沙漠 化 直接 因素 
与 间接 因素 比重 ,分 析 土 地 沙漠 化 的 驱动 缘由 。 


4 结果 与 分 析 


4.1 土地 沙漠 化 时 间 变 化 

对 2000—2018 年 塔里木 盆地 每 年 各 级 土地 沙 
漠 化 面积 进行 统计 ,如 图 2 所 示 。 由 图 可 知 ,2000 一 
2018 年 ,塔里木 盆地 植被 盖 度 呈 波 动 上 升 趋势 。19 a 


间 盆 地 内 植被 盖 度 年 栅 格 均值 增长 超过 50%。 非 
沙漠 化 面积 大 幅 增 长 ,累计 增长 20000 km ,增长 
1.67 倍 ,其 中 2005 2010, 2012 .2014 年 增长 较 快 ， 
2012 年 增长 最 快 ,增长 面积 达 9332.94 km’, JEWIN 
化 快速 增长 可 能 是 塔里木 河流 域 生态 输 水 量 增加 
和 退耕 还 林 还 草 政策 共同 作用 的 结果 呈 汪 。2006、 
2009 .2011.2015.2018 年 非 沙 漠 化 面积 则 出 现 小 幅 
下 降 , 最 大 下 降 面积 达 3000 km。 极 重度 沙漠 化 面 
积 整体 大 幅 下 降 ,累积 减少 20000 km*。2002 .2007、 
2012 .2015 年 极 重度 沙漠 化 面积 出 现 大 幅度 减少 ， 
其 中 2012 年 面积 减少 最 大 , 达 5384 km’, 而 2001、 
2004 2006 2016 年 则 出 现 小 幅 增 加 。 轻 度 沙漠 化 、 
中 度 沙 漠 化 .重度 沙漠 化 面积 变化 较 小 ,基本 呈 平 
4.2 土地 沙漠 化 空间 变化 
以 2000 年 为 基准 年 ,对 塔里木 盆地 进行 栅 格 累 
加 计算 ,结果 如 图 3 所 示 。 塔 里 木 盆地 土地 沙漠 化 整 
体 呈 改善 趋势 。 盆 地 改善 速率 为 44749.80 kma, 
恶化 速率 为 3000.52 km2.a-, 改 善 速率 远大 于 恶化 
速率 ,盆地 土地 沙漠 化 整体 呈 改 善 状态 。 改 善 主要 
以 低速 改善 和 中 速 改善 为 主 ,分别 为 20914.3 km?+ a’ 
和 19196.3 km .a ;其 次 为 高 速 改 善 和 极 高 速 改 善 ， 
分 别 为 3560.6 km-a” 和 1024.8 km…a'; 超 高 速 改 善 
较 小 , 仅 53.8 km …a'。 改 善 区域 广 泛 分 布 于 绿洲 内 
部 及 边缘 ,塔里木 河 下 游 沙 漠 化 改善 也 较为 明显 。 
恶化 速率 相对 较 低 ,主要 为 低速 恶化 ,速率 2265.2 
km2.a 5 其 次 是 中 速 恶化 ,633.6 km .a ;高 速 恶化 、 
极 高 速 恶化 . 超 高 速 恶 化 速率 较 小 。 恶 化 区 域 集中 
于 塔里木 河中 游 和 车 尔 臣 河 下 游 ,其 中 塔里木 河中 
游 ( 阿 拉 尔 至 恰 拉 水 库 段 ) 恶 化 最 为 严重 ,盆地 80% 
的 低速 .中 速 高速. 极 高 速 恶 化 均 集 中 于 此 ,恶化 


表 2 土地 沙漠 化 栅 格 累加 速率 分 级 


Tab.2 Classification of desertification rate used in the grid accumulation method 


分 级 超 高 速 恶化 极 高 速 恶化 高 速 恶 化 中 速 恶 化 低速 恶化 无 变化 
值 域 范 围 [-100.0 ,-60.0) [-60.0,-40.0) [-40.0 ,-26.7) [-26.7,-13.3) [-13.3,-6.0) [-6.0,6.0] 
分 级 低速 改善 中 速 改 善 高 速 改善 极 高 速 改 善 超 高 速 改善 

值 域 范围 (6.0,13.3] (13.3 ,26.7] (26.7 ,40.0] (40.0 ,60.0] (60.0, 100.0] 
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图 2 2000 一 2018 年 塔里木 盆地 土地 沙漠 化 变化 
Fig. 2 Desertification change in Tarim Basin during 2000—2018 


速率 达 2376 kwao ÆR EI PR CL BON HH 
显 ,主要 为 低速 .中 速 恶化 。 另 外 , 叶 尔 羌 河 与 喀什 
噶 尔 河 交 汇 处 出 现 超 高 速 恶化 区 ,不 容 忽 视 。 
4.3 土地 沙漠 化 转移 变化 

2000—2018 年 各 类 土地 沙漠 化 转移 面积 统计 
如 图 4 所 示 。 从 图 可 知 ,盆地 改善 转移 量 明显 高 于 
恶化 转移 量 。 改 善 转移 量 达 12179 km*…a ,恶化 转 
移 量 仅 8916 km .a'。 塔 里 木 贫 地 整体 呈 改 善 趋 
势 。 在 改善 转移 中 , 极 重 度 至 重度 转移 量 最 高 ,年 
均 转移 量 达 3920 km , 占 改 善 量 的 32.18%; 其 次 为 
轻 度 至 非 沙 漠 化 ,年 均 转 移 量 3290 km’, 占 改善 量 的 
27.01% . HERE 38 A BE .中 度 至 轻 度 转 移 量 相对 较 
小 ,转移 量 在 2400 km .a-: 左 右 , 分 别 占 改善 量 的 


21.18% 19.62%。 在 恶化 转移 中 ,重度 沙漠 化 至 极 
重度 沙漠 化 转移 量 最 大 ,年 均 转 移 量 达 3102 km’, 占 
总 恶化 面积 的 34.79%; 其 次 为 非 沙漠 化 至 轻 度 沙漠 
化 , 占 总 恶化 量 的 25.79%; 再 次 为 中 度 至 重度 ,转移 
量 2400 km …a'; 轻 度 至 中 度 , 转 移 量 1606 km2.a '; 
轻 度 至 重度 ,转移 量 仅 260 km? a Ati, EREK 
量 小 于 50 km…a ,未 在 图 4 中 表示 。 

各 转移 量 均 呈 年 际 波动 变化 。 波 动量 最 大 是 
轻 度 沙漠 化 到 非 沙 漠 化 ,振幅 达 8000 km ,其 次 是 极 
重度 至 重度 ,振幅 近 4000 km’; 再 次 是 非 沙 漠 化 至 轻 
度 .重度 至 极 重 度 ; 中 度 与 重度 、 轻 度 与 中 度 之 间 相 
互 转换 振幅 也 较 大 。 年 际 变 化 过 程 中 波峰 波 谷 相 
间 出 现 。 在 改善 转移 中 ,2003、2005、2007、2010、 
2012、2014、2017 年 出 现 波 峰 ,2012 年 达到 顶峰 ，; 
2004、2006、2009、2011、2013、2015、2018 年 出 现 波 
谷 ,波峰 波 谷 呈 相 间 分 布 。 恶 化 转移 与 改善 转移 相 
似 , 也 呈 波 动 状 态 ,但 恶化 转移 的 波峰 波 谷 与 改善 
转移 的 波 谷 波峰 几乎 相反 。 改 善 转移 的 波峰 对 应 
恶化 转移 的 波 谷 ,改善 转移 的 波 谷 则 对 应 于 恶化 转 
移 的 波峰 。 但 改善 转移 量 大 于 恶化 转移 量 ,盆地 的 
土地 沙漠 化 向 改善 方向 发 展 。 这 种 波动 变化 与 盆 
地 河流 来 水 .生态 输 水 存在 较 大 关系 1。 

为 直观 展示 各 转移 类 型 的 空间 分 布 ,选取 典型 
4 期 (2001、2006、2008、2010 年 ) 转 移 分 布 图 进行 展 
示 ( 图 5)。 从 转移 的 空间 分 布 上 看 , 轻 度 至 非 沙漠 
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图 3 2000 一 2018 年 塔里木 盆地 了 


Fig. 3 Accumulate spatial distribution of desertification using the grid accumulation method in Tarim Basin during 2000 一 2018 
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图 4 2000 一 2018 年 塔里木 盆地 土地 沙漠 化 转移 量 


Fig. 4 Transformed area of desertification in Tarim Basin during 2000—2018 
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图 5 塔里木 盆地 土地 沙漠 化 转移 空间 分 布 


Fig. 5 Spatial distribution of the transformed desertification in Tarim Basin 
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绿洲 区 内 ;中 度 至 轻 度 、 轻 度 至 中 度 ,重度 至 中 度 、 
中 度 至 重度 , 极 重 度 至 重度 .重度 至 极 重 分 布 相对 
集中 ,主要 分 布 于 塔里木 河中 上 游 . 车 尔 臣 河 下 游 、 
渭 干 河 下 游 . 叶 尔 芜 河 与 喀什 噶 尔 河 交汇 处 。 


5 讨论 


5.1 土地 沙漠 化 发 展 方向 

上 述 分 析 可 知 ,土地 沙漠 化 各 等 级 之 间 相 互 转 
换 , 且 转换 频繁 ,为 进一步 分 析 塔 里 木 贫 地 土地 沙 
漠 化 转移 方向 ,本 文 将 每 年 土地 沙漠 化 各 等 级 的 转 
移 来 源 统 计 如 图 6 所 示 。 非 沙漠 化 与 极 重 度 沙 漠 化 
的 转换 来 源 比 较 单 一 。 非 沙漠 化 94% 来 源 于 轻 度 
沙漠 化 , 极 重度 沙漠 化 96% 来 源 于 重度 沙漠 化 。 轻 
度 沙 漠 化 .中 度 沙 漠 化 .重度 沙漠 化 则 主要 由 2 种 十 
地 沙漠 化 转换 而 来 。 轻 度 沙漠 化 49% 来 自 中 度 沙 
漠 化 ,42% 来 自 非 沙漠 化 ;中 度 沙漠 化 57% 来 自重 度 
沙漠 化 ,36% 由 轻 度 沙漠 化 提供 ;重度 沙漠化 64% 由 
极 重 度 沙 漠 化 贡献 ,33% 由 中 度 沙 漠 化 贡献 。 从 转 
移 数量 上 看 ,塔里木 贫 地 的 土地 沙漠 化 转移 主要 发 
生 在 相 邻 的 土地 沙漠 化 等 级 之 间 ,转移 路 径 为 : 极 
重度 沙漠 化 一 重度 沙漠 化 一 中 度 沙漠 化 一 轻 度 沙 
漠 化 一 非 沙 漠 化 ,这 与 塔里木 盆地 土地 沙漠 化 改善 
现象 相 吻 合 。 

非 相 邻 的 土地 沙漠化 等 级 之 间 转 移 量 ( 跨 级 转 
移 量 ) 较 小 旦 相对 较 稳定 。 跨 级 转移 存在 2 个 明显 
阶段 。2001 一 2009 年 , 跨 级 转移 量 较 小 ,2010 一 
2018 年 转移 量变 大 。 在 跨 级 转移 中 , 跨 一 级 转移 量 
占 比较 多 ,如 极 重 度 与 中 度 、 重 度 与 轻 度 .中 度 与 非 
沙漠 之 间 的 相互 转换 ; 跨 2 级 及 2 级 以 上 转移 量 非 
常 少 。 可 见 塔 里 木 盆地 土地 沙漠 化 在 年 际 间 是 一 
个 逐步 连续 .逐年 波动 缓慢 变化 的 过 程 。 
5.2 土地 沙漠 化 成 因 

为 分 析 塔 里 木 盆 地 土地 沙漠 化 成 因 , 现 将 4.2 
部 分 中 的 栅 格 累加 计算 结果 按 土地 利用 类 型 (2010 
年 土地 利用 图 ) 进 行 统计 ,统计 结果 如 图 7 所 示 。 从 
统计 结果 可 知 , 土 地 沙漠 化 变化 主要 发 生 在 水 域 、 
草地 、 林 地、 耕地 和 未 利用 地 上 。 

土地 沙漠 化 恶化 区 主要 发 生 在 草地 .水域 、 林 
地 中 。 超 高 速 恶化 全 部 在 水 域 地 类 上 ; 极 高 速 恶 化 
主要 在 草地 与 水 域 地 类 上 ;高 速 和 中 速 恶化 主要 在 
草地 、 林 地 上 ;低速 恶化 则 主要 在 草地 、 建 设 用 地 、 


耕地 和 未 利用 地 上 。 从 地 类 上 看 ,草地 退化 是 沙漠 
化 恶化 的 主要 土地 类 型 ,其 次 是 林地 和 水 域 , 再 次 
是 耕地 建设 用 地 和 未 利用 地 。 

塔里木 盆地 草地 主要 分 布 于 河流 2 岸 和 尾 疗 湖 
周边 的 湿地 ,群落 种 类 单一 , 呈 隐 域 性 .中 生性 和 不 
稳定 性 59 , 极 易 受 河 水 和 地 下 水 影响 。 当 河水 和 地 
下 水 位 下 降 ,草地 迅速 响应 而 衰退 , 沙 地 裸露 ,沙漠 
化 加 重 汪 。 河 水 与 地 下 水 位 具有 年 际 波动 特性 ,从 
而 导致 草地 具有 年 际 波动 特性 ,这 可 能 是 塔里木 盆 
地 沙漠 化 年 际 波 动 的 主要 原因 。 塔 里 木 盆 地 林地 
主要 为 河岸 林 和 绿洲 防护 林 。 从 4.2.4.3 部 分 和 野 
外 验证 调查 可 知 ,林地 退化 区 域 主要 为 河岸 林 , 特 
是 远离 河道 的 天 然 林 , 退 化 主要 为 灌木 衰退 。 林 
地 对 河水 和 地 下 水 水 位 的 响应 速度 相 比 草地 要 慢 
一 些 , 因 此 林地 退化 主要 集中 在 中 速 恶 化 上 。 塔 里 
木 盆地 水 域 面积 受 人 类 活动 和 气候 变化 而 频繁 波 
动 ,变化 速度 迅速 ,因此 退化 速度 表现 为 超 高 速 。 
虽然 以 植被 盖 度 角度 分 析 得 出 水 域 退 化 不 一 定 能 
真实 反映 在 水 域 地 类 上 的 沙漠 化 变化 ,但 至 少 可 以 
反映 水 域 变 化 的 速度 。 塔 里 木 盆地 耕地 退化 也 不 
容 忽 视 ,在 高 速 、 中 速 和 低速 恶化 过 程 中 均 有 12% 
来 自 耕 田 , 耕 地 退化 主要 是 耕地 盐 渍 化 造成 的 耕地 
退化 和 弃 耕 形成 ,分 布 比较 分 散 。 另 外 ,在 低速 恶 
化 中 建设 用 地 占 近 20% ,建设 用 地 的 退化 原因 还 需 
要 进一步 深入 探寻 。 毕 上 所 述 ,塔里木 盆地 沙漠 化 
恶化 的 直接 原因 是 因 水 资源 变动 而 引起 交 ,特别 是 
河水 与 地 下 水 水 位 变动 造成 ,因此 水 是 塔里木 盆地 
沙漠 化 治理 的 关键 因子 1。 

土地 沙漠 化 改善 区 主要 发 生 在 林地 、 未 利用 地 
和 耕地 中 。 超 高 速 改 善 主要 发 生 在 未 利用 地 ; 极 高 
速 和 高 速 改善 主要 发 生 在 林地 ;中 速 和 低速 改善 则 
主要 发 生 在 林地 和 耕地 。 从 发 生地 类 可 知 , 林 地 是 
改善 的 主要 地 类 ;其 次 是 耕地 和 未 利用 地 ;再 次 是 
草地 。 与 土地 沙漠 恶化 相 比 ,草地 在 改善 过 程 中 占 
比 很 小 , 仅 占 极 高 速 和 高 速 改 善 10% 左 右 。 从 沙漠 
化 意识 调查 可 知 ,大 部 分 人 认为 治理 土地 沙漠 化 主 
要 途径 是 植树 造林 ”5。 塔 里 木 盆 地 也 一 样 ,在 绿洲 
内 部 和 边缘 进行 大 量 的 植树 造林 ,林地 得 到 改善 ; 
另外 , 塔 河 下 游 河 岸 林 因 塔 河 下 游 输 水 工程 地 持续 
实施 恢复 效果 也 较为 显著 ”。 绿 洲 防 护林 使 得 绿 
洲 内 部 生态 环境 改善 ,但 河岸 和 绿洲 外 围 的 荒漠 林 
却 依 然 在 退化 ,天 然 林 的 生态 环境 在 持续 恶化 。 塔 
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图 6 土地 沙漠 化 转化 来 源 


Fig.6 Sources of desertification transformation in Tarim Basin 
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里 木 盆 地 耕地 改善 占 比 也 较 大 。2002 年 
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图 7 各 速率 的 地 类 分 析 
Fig.7 Analysis of the desertification change rate of 


land-use type in Tarim Basin 


的 过 程 , 因 此 改善 速度 多 以 低速 和 中 速 改善 为 主 。 
未 用 地 改善 也 十 分 显著 ,在 极 高 速 改 善 中 占 比 较 
大 。 在 大 量 资金 支持 下 ,塔里木 盆地 开展 了 大 量 的 
治 沙 工程 ,如 且 末 76.7 km 河东 治 沙 工程 ,31.28 km? 
的 沙漠 公路 防风 林 工 程 。 这 些 治 沙 工程 使 未 利用 
地 植被 盖 度 迅速 提高 ,因此 未 用 地 改善 主要 为 极 高 
速 改 善 。 从 改善 地 类 和 改善 成 因 可 知 ,塔里木 盆 沙 
漠 化 改善 主要 为 人 为 因素 ,在 大 量 的 资金 投入 下 沙 
漠 化 得 以 大 面积 改善 。 这 种 改善 属于 外 来 输入 性 
改善 ,可 持续 性 较 差 ,一旦 资源 输入 停止 ,改善 区 域 
面临 恶化 而 恢复 到 原来 的 状况 。 塔 里 木 盆地 沙漠 
化 治理 的 手段 和 方式 还 需要 进一步 深入 研究 。 
5.3 栅 格 累加 分 析 法 

从 4.1 和 4.3 部 分 可 知 , 土 地 沙漠 化 是 一 个 长 时 
间 波 动 变化 过 程 ,这 种 波动 性 可 能 是 土地 沙漠 化 进 
程 中 的 一 种 常态 , 若 不 考虑 这 种 波动 性 仅 选取 几 期 
遥感 影像 进行 分 析 ,其 结果 往往 具有 偶然 性 ,导致 
结果 可 能 存在 偏差 ,甚至 得 出 完全 相反 的 结论 , 同 
时 波动 变化 过 程 中 往往 存在 一 些 的 累积 效应 所 ,这 
种 累积 效应 通常 预示 了 一 定 的 沙漠 化 发 展 方向 , 累 
耻 效 应 分 析 结 果 可 为 沙漠 化 演变 研究 提供 帮助 。 
因此 ,合理 的 土地 沙漠 化 变化 监测 需要 考虑 沙漠 化 
的 波动 性 与 累积 效应 。 本 文 使 用 栅 格 累加 法 利用 
多 年 连续 的 遥感 数据 进行 栅 格 累积 分 析 , 既 去 除 波 
动产 生 的 误差 又 分 析 了 累加 状态 ,以 此 得 到 较为 真 
实 的 生态 环境 演变 方向 演变 趋势 和 变化 速率 。 如 


为 生态 环境 监测 提供 了 一 种 有 效 地 监测 方法 与 

本 文 基 于 栅 格 累加 法 分 析 塔 里 木 贫 地 的 土地 
沙漠 化 变化 趋势 和 变化 空间 ,定量 确定 了 塔里木 盆 
地 土地 沙漠 化 转移 空间 分 布 、 转 移 过 程 和 转移 方 
向 ,对 土地 沙漠 化 变化 有 了 更 直观 、 更 深刻 的 理解 ， 
HAT aR SAAS ESTA, HA AS: 
合 土地 利用 类 型 定量 给 出 塔里木 盆地 土地 沙漠 化 
的 原因 ,为 沙漠 化 治理 提供 了 准确 的 空间 和 成 因 
信息 。 

但 本 文 仅 考 虑 了 NDVI 这 一 指标 ,对 气候 E 
壤 、 人 类 活动 等 探讨 还 远 远 不 够 ,还 需 进 一 步 改 
进 。 另 外 ,本文 使 用 的 数据 分 辩 率 过 低 ,更 高 分 辩 
率 数据 也 值得 探讨 。 


6 结论 


本 文 以 塔里木 盆地 2000—2018 年 MODIS-ND- 
VI 1 km 月 植被 数据 为 基础 ,运用 栅 格 累加 法 分 析 塔 
里 木 盆地 土地 沙漠 化 变化 ,得 出 以 下 结论 : 

(1) 2000—2018 年 塔里木 盆地 土地 沙漠 化 是 一 
个 年 际 波 动 变化 过 程 。 盆 地 内 各 等 级 土地 沙漠 化 
之 间 相 互 转 换 频繁 ,年 际 转换 面积 呈 显 著 波动 , 造 
成 土地 沙漠 化 面积 年 际 波动 明显 。 

(2) 栅 格 累加 法 原理 简单 过程 明了 ,能 克服 年 
际 波动 性 产生 的 误差 ,为 连续 长 时 间 空 间 数 据 的 演 
变 趋势 演变 速率 和 空间 分 布 提 供 了 一 种 简单 实用 
的 新 方法 。2000 一 2018 年 塔里木 盆地 土地 沙漠 化 
分 析 揭 示 :盆地 土地 沙漠 化 整体 呈 改 善 趋势 ;改善 
速率 44749.80 km2.a ,改善 以 低速 改善 和 中 速 改善 
为 主 ,改善 区 域 主要 分 布 在 绿洲 内 部 和 边缘 ;恶化 
速率 3000.52 km,a ,主要 为 低速 恶化 ;恶化 区 域 主 
要 分 布 于 塔里木 河中 游 和 车 尔 臣 河 下 游 。 

(3) 栅 格 累加 法 结合 转移 和 矩阵 分 析 为 定量 获取 
土地 沙漠 化 转移 路 径 转移 数量 转移 空间 分 布 提 
供 了 一 种 新 的 手段 。 本 文 首次 定量 地 获取 塔里木 
盆地 土地 沙漠 化 的 转移 路 径 ,其 转移 路 径 为 : 极 重 
度 沙 漠 化 一 重度 沙漠 化 一 中 度 沙漠 化 一 轻 度 沙漠 
化 一 非 沙漠 化 ;其 中 轻 度 沙 漠 化 与 非 沙 漠 化 之 间 的 
转移 量 最 大 ,主要 发 生 于 绿洲 区 内 ,空间 分 布 较 广 ; 
中 度 与 轻 度 、 重 度 与 中 度 、 极 重度 与 重度 之 间 的 转 
换 量 相 对 较 小 ,分 布 相对 集中 ,主要 分 布 于 塔里木 
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(4) 土地 利用 类 型 结合 栅 格 累加 法 为 土地 沙漠 


成 因 定量 分 析 提 供 了 一 种 新 的 和 尝试。 本文 分 析 显 
示 塔 里 木 岔 地 土地 退化 主要 是 草地 退化 ,改善 主要 
为 林地 和 耕地 改善 。 林 地 以 绿洲 内 部 和 外 部 防风 
林 改 善 为 主 , 耕 地 改善 则 主要 是 盐 渍 化 改造 和 种 植 
结构 调整 形成 。 改 善 主 要 是 因 人 类 主动 活动 形成 。 


致谢 :感谢 https://search.earthdata.nasa.gov/ 网 站 


提供 MODIS NDVI 数 据 ,http://www.resdc.cn/ 网 站 提 


供 土 地 利用 数据 。 
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刘 WAS: 基于 栅 格 累加 法 的 塔里木 盆地 沙漠 化 分 析 


Desertification analysis based on grid accumulation method in 


Tarim Basin, China 


LIU Chang”，LI Chengzhi'’’, LI Shenghui'’, FU Xiaolei”, SHI Qingdong 


(1. Institute of Arid Ecology & Environment of Xinjiang University, Urumqi 830046, Xinjiang, China; 2. College of Resources 
and Environmental Sciences, Xinjiang University, Urumqi 830046, Xinjiang, China; 3. Key Laboratory of Oasis Ecology 
Ministry of Education, Xinjiang University, Urumqi 830046, Xinjiang, China) 


Abstract: Currently, sandy desertification is one of the most serious environmental problems worldwide. It is 
also the most crucial ecological environmental problem in the arid areas in western China. The Tarim Basin, 
located in northwest China, is extremely arid, with little precipitation, high evaporation, and sparse vegetation. 
This basin includes the second largest flowing desert, and the problem of sandy desertification is extremely 
serious. In this paper, the desertification was monitored yearly using the data that were synthesized from the 
monthly data of MODIS-NDVI 1 km from 2000 to 2018 using the maximum synthesis method in the Tarim Basin. 
On the basis of the grid accumulation method, change rate, change trend, and spatial distribution, the causes of 
the sandy desertification were analyzed in the Tarim Basin. The results showed that: (1) the process of sandy 
desertification from 2000 to 2018 in the Tarim Basin fluctuated annually, and if the deviation calculation method 
had been used to monitor the desertification with several remote sensing images, it would have obtained the 
wrong results. (2) The grid accumulation method can quantitatively analyze the change rate, evolution trend, and 
spatial characteristics of sandy desertification, and the method is suitable for the analysis of the process of the 
interannual fluctuation. (3) The new method, combining the grid accumulation method with the transfer matrix 
analysis method, was presented. This method was used to quantitatively obtain the transfer path, transfer quantity, 
and transfer spatial distribution of the sandy desertification in the Tarim Basin. It revealed that the transfer path of 
the sandy desertification in the Tarim Basin changed from extremely severe desertification to severe 
desertification, to moderate desertification, to mild desertification, and to non-desertification. Generally, the 
desertification trend in the Tarim Basin was from extremely severe desertification to non- desertification. The 
improved areas were distributed in the oases, but the deteriorated areas were mainly located in the middle reaches 
of the Tarim River and lower reaches of the Qarqan River. (4) The causes of the desertification were 
quantitatively analyzed using the method that combined the land-use types and grid accumulation method. This 
method found that the most important cause of desertification was human activities, the main type of degradation 
was the degradation of the grassland, and the main types of improvement were those of forest land and cultivated 
land. This research provided a simple and easy method for the dynamic monitoring of desertification. Its results 
can accurately show the change rate, evolution trend, spatial distribution, and causes of land desertification, of 
which accurate information is necessary to control the desertification of the land. 
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